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Summary 

As  a  direct  response  to  the  need  for  further  performance 
gains  from  current  multistage  axial  compressors,  an 
investigation  of  advanced  aerodynamic  design  concepts 
that  will  lead  to  compact,  high-efficiency,  and  wide- 
operability  configurations  is  being  pursued.  Part  I  of  this 
report  describes  the  projected  level  of  technical  advance¬ 
ment  relative  to  the  state  of  the  art  and  quantifies  it  in 
terms  of  basic  aerodynamic  technology  elements  of 
current  design  systems.  A  rational  enhancement  of  these 
elements  is  shown  to  lead  to  a  substantial  expansion  of  the 
design  and  operability  space.  Aerodynamic  design  con¬ 
siderations  for  a  four-stage  core  compressor  intended  to 
serve  as  a  vehicle  to  develop,  integrate,  and  demonstrate 
aerotechnology  advancements  are  discussed.  This  design 
is  biased  toward  high  efficiency  at  high  loading.  Three- 
dimensional  blading  and  spanwise  tailoring  of  vector 
diagrams  guided  by  computational  fluid  dynamics  (CFD) 
are  used  to  manage  the  aerodynamics  of  the  high-loaded 
endwall  regions.  Certain  deleterious  flow  features,  such 
as  leakage-vortex-dominated  endwall  flow  and  strong 
shock-boundary-layer  interactions,  were  identified  and 
targeted  for  improvement.  However,  the  preliminary 
results  were  encouraging  and  the  front  two  stages  were 


extracted  forfurther  aerodynamic  trimming  using  a  three- 
dimensional  inverse  design  method  described  in  part  II 
of  this  report. 

The  benefits  of  the  inverse  design  method  are  illus¬ 
trated  by  developing  an  appropriate  pressure-loading 
strategy  for  transonic  blading  and  applying  it  to  reblade 
the  rotors  in  the  front  two  stages  of  the  four-stage  configu¬ 
ration.  Multistage  CFD  simulations  based  on  the  average 
passage  formulation  indicated  an  overall  efficiency 
potential  far  exceeding  current  practice  for  the  front  two 
stages.  Results  of  the  CFD  simulation  at  the  aerodynamic 
design  point  are  interrogated  to  identify  areas  requiring 
additional  development.  In  spite  of  the  significantly  higher 
aerodynamic  loadings,  advanced  CFD-based  tools  were 
able  to  effectively  guide  the  design  of  a  very  efficient 
axial  compressor  under  state-of-the-art  aeromechanical 
constraints. 


Introduction 

The  intensely  competitive  aeropropulsion  sector  places 
some  of  the  most  stringent  requirements  on  turbo¬ 
compressors.  A  common  challenge  is  aerothermo- 
dynamic  optimization  subject  to  conflicting  economic 
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and  system  constraints,  such  as  reduced  cost,  low  environ¬ 
mental  emissions,  and  wide  operability  under  adverse 
conditions.  An  interesting  account  of  the  first  50  years  of 
aeropropulsion  gas  turbines  was  offered  by  Singh  (ref.  1 ), 
who  clarified  the  trail  taken  to  improve  fuel  economy  by 
way  of  higher  thermal  and  propulsive  efficiencies. 
Figure  1  shows  the  effects  of  component  efficiency, 
compressor  delivery  temperature,  and  turbine  entry  tem¬ 
perature  on  the  thermal  efficiency  of  a  simple  Brayton 
cycle  using  a  hydrocarbon  fuel  and  metallic  materials. 
The  implication  is  that  overall  polytropic  efficiencies  of 
94  percent  with  compression  ratios  in  the  range  of  60: 1 
and  beyond  must  be  realized  to  arrive  at  thermal  efficien¬ 
cies  surpassing  55  percent.  Techniques  for  achieving 
low  losses  at  high  aerodynamic  loadings  have  therefore 
received  renewed  interest  in  many  recent  turbomachinery 
studies  (refs.  2  and  3). 

Turbocompression  technology  has  been  advanced 
continuously  by  higher  work  capacity  per  stage  as  a  result 
of  increases  in  rotor  speed,  aerodynamic  loading,  and 
throughflow  Mach  numbers.  Using  sophisticated  diag¬ 
nostic  tools  involving  CFD  and  measurement  tech¬ 
niques,  more  suitable  blade  shapes  with  relatively  low 
losses  at  higher  diffusion  and  Mach  number  levels  have 
been  deployed.  Better  materials  and  matured  structural 
analysis  methods  have  also  allowed  increases  in  rotor 
speed  with  significant  weight  reductions.  Improved 
mechanical  design  and  fabrication  techniques  have 
raised  the  quality  of  current  products. 


The  quest  for  further  aerodynamic  performance 
advancements  is  becoming  progressively  more  difficult 
because  of  a  dwindling  residue  of  losses.  Prohibitive 
demands  such  as  performance  invariance  to  various  oper¬ 
ating  conditions  and  hardware  degradation  are  imposing 
severe  aerodynamic  limitations.  Easing  these  limitations 
requires  the  application  of  new  technology  to  manage  the 
particular  flow  structures  responsible  for  performance 
shortfalls.  Effective  strategies  for  achieving  this  are  not  at 
all  clear.  A  vital  issue  is  whether  to  advance  by  refining 
and  extending  well-proven  concepts,  perhaps  in  the  face 
of  diminishing  returns,  or  by  changing  to  something 
more  unpredictable  yet  inviting. 

A  research  program  being  conducted  at  the  Glenn 
Research  Center  is  investigating  advanced  design  con¬ 
cepts  that  will  lead  to  compact,  high-efficiency,  and  wide- 
operability  compressors.  This  is  a  direct  response  to  the 
need  for  further  performance  gains  from  current 
turbomachinery  systems.  To  service  these  gains,  a  com¬ 
bination  of  evolutionary  and  revolutionary  approaches  to 
technology  development  was  selected.  The  evolutionary 
approach  employs  advancements  in  simulation  techniques 
to  refine  traditional  design  concepts  in  a  bid  for  higher 
efficiencies  at  increased  aerodynamic  loading  levels, 
whereas  the  revolutionary  approach  attempts  to  explore 
unconventional  concepts  and  paradigms  for  increased 
pressure  ratio,  higher  efficiencies,  and  wider  operability. 


Turbine  entry 


Figure  1. — Impact  of  component  efficiency  on  thermal  efficiency  (Singh,  ref.  1) 
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Part  I  of  this  report  depicts  the  level  of  technical 
advancement  being  sought,  quantifies  this  advancement 
in  terms  of  the  basic  aerodynamic  technology  elements  of 
current  design  systems,  and  identifies  a  four-stage  core 
compressor  configuration  for  advanced  technology 
development.  A  preliminary  aerodynamic  design  of  the 
four-stage  configuration  is  described  in  the  context  of 
current  practice.  Custom-tailored  blade  contours  were 
initially  developed  from  mean  camber  surfaces  derived 
from  axisymmetric  stream  surfaces  generated  by  an 
inverse  throughflow  calculation.  Isolated  blade-row 
CFD  simulations  provided  some  guidance  to  fine-tune 
the  blading.  A  multistage  CFD  computation  using  a 
mixing  plane  scheme  indicated  an  efficiency  potential  for 
the  four-stage  compressor  that  is  within  current  stan¬ 
dards.  In  addition,  performance-limiting  flow  features, 
such  as  leakage-vortex-dominated  endwall  action  and/or 
reaction  and  strong  shock-boundary-layer  interactions, 
were  identified  and  targeted  for  improvement.  The  possi¬ 
bility  of  further  performance  gains  requires  a  design 
method  that  effectively  extends  the  custom-tailored 
blading  philosophy  into  the  endwall  regions  and 
adequately  allows  matching  of  blade  rows  in  the  multi¬ 
stage  environment. 

In  part  II  of  the  report,  an  aerodynamic  design  of  the 
front  two  stages  of  the  four-stage  configuration  is  devel¬ 
oped  using  the  best  available  CFD  tools.  The  three- 
dimensional  inverse  design  method  of  Dang  (ref.  4)  was 
used  to  define  advanced  transonic  rotor  blading.  To 
illustrate  the  benefits  of  this  method,  the  front  two  stages 
of  the  four-stage  configuration  were  rebladed  and  the 
performance  potential  evaluated.  The  multistage  aerody¬ 
namic  environment  was  modeled  using  the  substantiated 
method  of  Adamczyk  (ref.  5).  This  modeling  permitted  a 
reliable  prediction  of  the  performance  potential  of  the 
two-stage  group  and  guided  design  revisions,  including 
the  facilitation  of  stage  matching.  The  results  show  the 
critical  role  of  this  method  and  illustrate  achievable 
performance  using  an  evolutionary  approach  buttressed 
by  advanced  turbomachinery  CFD. 

Geometric  information  for  all  blade  rows  is  given  in 
appendix  A,  and  span  wise  profiles  of  axisymmetric  aver¬ 
aged  flow  properties  extracted  from  the  multistage  CFD 
simulation  of  the  front  two-stages  at  the  aerodynamic 
design  point  are  given  in  appendix  B. 


Symbols 

Aex  passage  exit  area,  in. 

Aq  tangential  projection  of  blade  surface  area 

Dfador  diffusion  factor 

E  effectivity  ratio 

/  blade  mean  camber  surface,  f(r,z) 

/stack  blade  stacking  line>/stack(r) 

AH  total  enthalpy  change,  BTU/lbm 

Ldiff  equivalent  diffusion  length 

Ldlff/^An  ratio  of  equivalent  diffusion  length 
to  effective  passage  exit  width 
M|  inlet  Mach  number 

rir  mass  flow  rate,  lbm/s 

P  total  pressure,  psia 

PR  pressure  ratio 

p  static  pressure,  psia 

A p  static-pressure  difference  between  blade 

upper  and  lower  surfaces  at  fixed  axial 
positions,  A p  =  A p(r,z) 

Rec  chord  Reynolds  number 

T  blade  thickness  distribution,  T(r,z) 

tig  ratio  of  clearance  to  staggered  spacing 

U  blade  speed,  ft/s 

V  velocity,  ft/s 

Vq  absolute  tangential  velocity,  ft/s 

Wc  inlet  corrected  flow,  lbm/s 

z/s  ratio  of  axial  spacing  to  pitch 

a  swirl  angle  from  axial  direction,  deg 

r|  efficiency,  percent 
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A 

blade  aspect  ratio;  ratio  of  average  length  to 

pitchline  chord 

a 

solidity;  ratio  of  chord  to  average  tangential 

spacing 

<P 

flow  coefficient,  V_  !Um 

V 

work  coefficient,  A  HI  u}n 

Q 

corrected  mechanical  speed,  rpm 

Subscripts: 

abs 

absolute 

ex 

exit 

in 

inlet 

le 

leading  edge 

m 

meanline 

max 

maximum 

P 

polytropic 

ref 

reference 

rel 

relative 

t 

total 

te 

trailing  edge 

tip 

tip 

tt 

total  to  total 

z 

axial 

Part  I:  Level  of  Technical  Advancement 

Basic  Aerodynamic  Technology  Elements 

Aerodynamic  loading  and  loss  synthesis, — Turbo¬ 
compressor  aerodynamics  comprises  a  complex  array  of 
intimately  linked  thermofluid  processes  that  are  difficult 
to  separate  into  neat  compartments.  However,  simplified 
models,  which  sufficiently  incorporate  the  important 
physical  phenomena,  have  been  synthesized  into  design 
systems.  Three  essential  aerodynamic  technology  ele¬ 
ments  of  such  design  systems  are  ( 1 )  entropy  production 
or  loss;  (2)  aerodynamic  transport  including  spanwise 


mixing;  and  (3)  aerodynamic  stability.  The  demarcation 
of  these  three  elements  is  sometimes  distorted  by  various 
empirical  correlations  that  implicitly  account  for  aerody¬ 
namic  transport.  Denton  (ref.  6)  gives  a  good  physical 
description  of  entropy  production  in  turbomachines  and 
Greitzer  (ref.  7)  does  likewise  for  aerodynamic  stability. 
Adamczyk  (ref.  5)  and  Adkins  and  Smith  (ref.  8)  discuss 
aerodynamic  transport  phenomena  in  turbocompressors. 

An  adaptation  of  the  loss  model  of  Koch  and  Smith 
(ref.  9)  and  the  stage  maximum  static-pressure-rise  pre¬ 
diction  method  of  Koch  (ref.  10)  was  used  to  estimate 
both  the  efficiency  potential  and  the  maximum  aerody¬ 
namic  loading  capability  for  the  present  study.  The  losses 
are  grouped  according  to  source:  (1)  viscous  dissipation 
along  the  blade-surface  profile  and  the  resultant  wake 
mixing;  (2)  endwall  viscous  shear  including  interactions 
with  secondary  and  leakage  flows;  and  (3)  shocks  on  the 
blading.  Blade-surface-profile  losses  are  linked  to 
suction-surface  diffusion  (V^ax/V^,),  blade  trailing-edge 
thickness,  Reynolds  number,  and  surface  roughness. 
Endwall  losses  that  have  been  determined  from  hub  and 
casing  viscous  layer  measurements  are  related  to  aspect 
ratio,  solidity,  blade  stagger,  endwall  clearance,  blade- 
row  axial  spacing,  and  aerodynamic  loading  level.  The 
shock  loss  model  relates  passage  shock  losses  to  inlet  and 
exit  Mach  numbers.  A  leading-edge  bow  shock  loss  is 
calculated  based  on  the  inlet  Mach  number  and  leading- 
edge  thickness. 

The  stage  maximum  static-pressure-rise  potential  is 
linked  to  stage  geometric  parameters  (e.g.,  solidity, 
aspect  ratio,  stagger,  etc.)  by  a  correlation  based  on  an 
analogy  between  compressor  blade  passages  and  straight 
diffusers.  There  are  also  corrections  for  Reynolds  num¬ 
ber,  endwall  leakage,  and  blade-row  axial  spacing.  The 
blading  is  averaged  over  a  stage  without  separating  the 
rotor  and  stator  performance.  This  is  a  simple  yet  effec¬ 
tive  method  to  complement  conceptual  design  calcula¬ 
tions  made  at  a  representative  section  of  a  machine  (i.e., 
meanline  with  implicit  accounting  of  aerotransport).  A 
relative  aerodynamic  loading  parameter  is  defined  by  the 
ratio  of  the  actual  stage  pressure-rise  coefficient  to  that 
predicted  by  the  correlation  (the  effectivity  ratio  E). 
Hence,  an  effectivity  of  1  (E=  1 .0)  indicates  limit  loading. 
It  can  be  argued  that  this  is  a  more  appropriate  indicator 
of  aerodynamic  loading  from  the  point  of  view  of  insta¬ 
bility  initiated  by  the  endwall  flow.  Measurements 
(ref.  1 1 )  have  indicated  that  for  well-designed  blading 
with  adequate  circulation  capacity,  aerodynamic  instabil¬ 
ity  is  likely  to  occur  when  endwall  viscous  layers  increase 
to  a  limiting  value  dependent  on  effectivity  and  clearance. 
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The  loss  and  aerodynamic  loading  models  are  unified  in 
the  sense  that  endwall  loss  and  aerodynamic  blockage  are 
interrelated  through  aerodynamic  loading.  Thus,  it  is 
possible  to  determine  the  maximum  pressure  rise  that  can 
be  obtained  with  moderate  losses  but  not  necessarily  to 
determine  the  point  of  instability. 

Anticipated  technology  advancements. — Changes  to 
the  loss  and  peak  pressure-rise  models  were  made  to 
reflect  anticipated  aerodynamic  technology  advancements. 
It  was  surmised  that  low  blade  profile  losses  could 
be  extended  to  higher  diffusion  levels  by  employing 
boundary-layer  control  techniques.  Cascade  tests  and 
boundary-layer  analysis  have  shown  the  potential  for 
effecting  this  through  diffusion  control  using  tailored 
contours,  staged  diffusion  (e.g.,  tandem  and/or  splittered 
architectures),  or  active  diffusion  control  (e.g.,  surface 
transpiration  or  surface  morphing).  Figure  2  shows  the 
postulated  improvement  in  the  ratio  of  the  blade  trailing- 
edge  momentum  thickness  to  the  chord  (i.e.,  viscous 
dissipation)  relative  to  the  state  of  the  practice  for  a  fixed- 
chord  Reynolds  number  Rec,  an  inlet  Mach  number  Mv 
and  an  axial  velocity-density  ratio  (AVDR).  Note  that 
the  reference  is  the  state  of  the  practice  extracted  from 
the  method  of  Koch  and  Smith  rather  than  state  of  the  art 
because  there  are  some  uncertainties  as  to  what  the 
latter  is.  Also  indicated  is  a  data  point  assembled  from 
measurements  (ref.  12)  on  a  compressor  stator  equipped 
with  bleed  slots  for  boundary-layer  control.  Research 


Projections  with  advanced 


Suction-surface  diffusion  ratio,  Vmax/Vfe 


Figure  2. — Enhancements  to  profile  loss  model.  Fixed- 
chord  Reynolds  number,  Rec,  106;  inlet  Mach  number. 
Mi,  0.05;  axial  velocity-density  ratio,  AVDR,  1.0. 


currently  underway  should  provide  the  knowledge 
required  for  practically  achieving  this  control. 

The  stage  maximum  static-pressure-rise  capability  was 
postulated  to  increase  relative  to  the  state  of  the  practice 
as  shown  in  figure  3  for  a  fixed  Reynolds  number  Re,  the 
ratio  of  clearance  to  staggered  spacing  tig,  and  the  ratio  of 
axial  spacing  to  pitch  zJs.  The  abscissa,  LdjffA/Aex  » is 
the  ratio  of  an  equivalent  diffusion  length  Ldiff  to  an 
effective  passage  exit  width  (square  root  of  passage  exit 
area  Aex).  This  ratio  can  be  expressed  in  terms  of  concep¬ 
tual  design  parameters  (solidity,  aspect  ratio,  stagger)  at 
a  representative  section  of  the  machine.  Shown  on  the 
ordinate  is  the  stage  maximum  effective  static-pressure- 
rise  coefficient  averaged  over  the  rotor  and  stator  (see 
ref.  10  for  details).  Also  indicated  in  the  figure  is  a  data 
point  from  a  low-speed  rig  with  high-stagger,  forward- 
swept  blading.  This  data  point  is  clearly  above  the  state- 
of-the-practice  level  and  reasons  for  this  are  under 
investigation  using  a  simple  physical  model  developed 
by  Khalid  et  al.  (ref.  13)  for  clearance-related  blockage. 
The  results  from  this  model  provide  a  quantification  of 
aerodynamic  blockage  and  a  framework  for  screening 
endwall  flow-management  strategies  based  on  casing 
treatment,  vector  diagrams,  or  pressure-loading  distribu¬ 
tion.  The  model  corroborates  the  experiments  of  Lee  and 
Greitzer(ref.  14),  who  demonstrated  significant  enhance¬ 
ments  to  the  peak  static-pressure  rise  of  a  blade  row  by  the 
downstream  removal  of  low-momentum  flow  and  the 


Equivalent  conical  diffusion  angle  parameter, 
fdiff/^ex)0-5 

Figure  3. — Enhancements  to  stage  maximum  static- 
pressure  rise.  Fixed  Reynolds  number,  Re,  1 30  000; 
ratio  of  clearance  to  staggered  spacing,  tig,  0.055; 
ratio  of  axial  spacing  to  pitch,  z/s,  0.38. 
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upstream  injection  of  flow  with  sufficient  streamwise 
momentum. 

It  was  judged  that  managing  the  endwall  flow  (e.g., 
secondary  and  leakage  flow  control)  would  also  allow  the 
blading  to  be  more  effective  in  that  region,  thereby 
reducing  the  required  torque  supplied  to  the  rotor.  This  is 
accounted  for  in  the  endwall  loss  model  by  a  17-percent 
increase  in  the  so-called  tangential  force  thickness  as 
described  by  Smith  (ref.  11).  All  other  loss  sources, 
including  direct  shock  losses,  remained  unmodified.  The 
loss  and  aerodynamic  loading  models  were  integrated  in 
a  computerized  preliminary  design  and  analysis  proce¬ 
dure,  which  allowed  perturbations  to  specific  aerody¬ 
namic  technology  elements  to  be  readily  appraised. 

Extending  the  Design  Envelope 

The  consequences  of  the  aerodynamic  advancements 
previously  discussed  were  evaluated  by  comparing  the 
design  space  of  a  prototypical  stage  with  and  without 
these  advancements.  A  so-called  Smith  chart  is  used  to 
represent  the  design  and  operational  space  of  a 
turbocompressor  stage  with  a  given  geometric  form  and 
flow  type.  This  space  is  spanned  by  a  suitably  averaged 
flow  coefficient  (cp  =  VJUW)  and  work  coefficient 
(\| /  =  AH/U 2).  A  relative  aerodynamic  loading  para¬ 
meter  E  (V|/,  cp,  geometric  form,  flow  type)  and  an 
efficiency  q  (v|/,  cp,  geometric  form,  flow  type)  can  be 
superimposed  on  this  chart.  It  should  be  noted  that  the 
work  coefficient  is  generally  not  a  direct  measure  of 
aerodynamic  loading  as  defined  herein  because  it  has  no 
explicit  link  to  mean  static-pressure  rise.  This  is  reflected 
in  the  additional  functional  dependence  of  loading  on 
flow  coefficient  and  geometric  form.  The  basic  design 
challenge  is  to  generate  geometry  that  establishes  a  flow 
satisfying  the  \|/-cp  requirements  at  high  efficiency  with¬ 
out  exceeding  the  relative  aerodynamic  loading  limit. 

Geometric  parameters  for  a  prototypical  single- 
stage  axial  compressor  with  a  21  -in.  maximum  diameter 
were  selected  to  illustrate  the  expanded  design  envelope 
and  enhanced  performance  trends.  A  Smith  chart  for 
these  parameters  is  shown  in  figure  4(a)  using  state-of- 
the-practice  aerodynamic  technology  and  in  figure  4(b) 
with  the  postulated  technology  advancements.  These 
plots  are  based  on  the  following  design  choices:  (1)  zero 
swirl  at  stage  inlet;  (2)  60-percent-rotor-area  contraction; 
(3)  constant  radius  casing;  and  (4)  axial  velocity  ratio  of 
1.0  across  the  stage.  Since  the  rotor  inlet  specific  flow 
(40.0  lbm/s-ft2)  is  fixed,  varying  the  flow  coefficient 


Effectivity 


Flow  coefficient,  <p  =  Vz/Um 

Figure  4. — Expanded  design  envelope  and  enhanced 
performance  trends  for  prototypical  single-stage  axial 
compressor  with  21  -in.  maximum  diameter.  Geometric 
parameters:  rotor  aspect  ratio,  0.940;  stator  aspect 
ratio,  1.250;  rotor  solidity,  1.800;  stator  solidity,  1.700; 
Reynolds  number,  3.26x1 06;  ratio  of  clearance  to 
staggered  spacing,  0.0122.  (a)  Aerodynamic  design 
space  for  state-of-the-practice  technology,  (b)  Ex¬ 
tended  design  envelope  with  advanced  technology. 


sweeps  through  both  subsonic  and  transonic  stages.  In 
addition,  fixing  the  inlet  swirl  and  rotor  area  contraction 
generates  a  range  of  stage  reactions  as  the  flow  coefficient 
is  varied.  The  design  point  for  individual  stages  of  this 
particular  type  may  be  placed  anywhere  below  the  bound- 
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Rotor  tip  speed,  Utip.  ft/s 


Effectivity 


Rotor  tip  speed,  l/tip,  ft/s 

Figure  5.— Performance  trends  for  prototypical  stage, 
(a)  State-of-the-practice  technology,  (b)  Postulated 
advanced  technology. 


ary  E  =  1.0.  For  a  given  relative  aerodynamic  loading 
level  E,  flow  and  work  coefficients  can  be  selected  to 
maximize  efficiency. 

Evident  from  figure  4  is  a  significant  extension  of  the 
design  space  in  the  direction  of  higher  work  coefficients 
made  possible  by  the  postulated  aerodynamic  technology 
advancements.  A  potential  increase  of  three  points  in 
peak  polytropic  efficiency  is  depicted  along  with  the 
availability  of  more  options  for  high-efficiency  designs. 
Although  the  operability  space,  including  flow  range,  is 
not  directly  addressed  herein,  the  choice  of  how  far  and  in 
which  direction  to  extend  the  design  space  is  greatly 
influenced  by  operability  considerations.  To  illustrate 


specific  performance  trends,  figure  5  presents  the  effect 
of  rotor  tip  speed  and  stage  aerodynamic  loading  on  stage 
total-pressure  ratio  and  polytropic  efficiency  for  the  pro¬ 
totypical  stage.  Figure  5(a)  is  based  on  state-of-the- 
practice  technology  whereas  figure  5(b)  used  the  postulated 
aerodynamic  technology  advancements.  The  curves  indi¬ 
cate  that  higher  stage  pressure  ratios  can  be  obtained  by 
either  increased  tip  speed  or  increased  stage  loading. 
Maximum  dividends  occur  when  both  high  loading  and 
high  tip  speed  are  selected.  The  attainment  of  gains  in 
pressure  ratio  is  complicated  by  the  need  for  high 
efficiency  and  wide  operability.  As  loading  is  increased, 
losses  tend  to  increase,  and  the  potential  for  flow  break¬ 
down  and  aerodynamic  instability  is  much  greater.  For 
the  range  of  blade  speeds  shown  in  figure  5,  the  rotor  tip 
relative  Mach  numbers  vary  from  subsonic  to  high  super¬ 
sonic  (i.e.,  greater  than  1 .6)  values  that  can  result  in  severe 
shock-related  viscous  losses. 

For  a  fixed  geometric  envelope,  the  anticipated  aerody¬ 
namic  technology  advancements  suggest  the  potential  for 
substantial  expansion  of  the  design  space  in  the  direction 
of  higher  work  coefficients  and  stage  total-pressure  ratio. 
In  addition,  high  efficiencies  are  extended  to  more  el¬ 
evated  loadings.  An  increased  pressure  ratio  per  stage  has 
traditionally  been  realized  partly  because  of  higher  rotor 
speeds  and  primarily  because  of  a  lower  aspect  ratio, 
increased  solidity,  and  higher  stagger  blading.  These 
design  choices,  as  can  be  inferred  from  figure  3,  are  often 
made  to  increase  stage  maximum  static-pressure  rise, 
thereby  insuring  adequate  aerodynamic  stability,  which 
can  often  be  accompanied  by  increased  losses.  A  report 
by  Wisler,  Koch,  and  Smith  (ref.  15)  illustrates  the 
interactions  between  different  design  choices  and  offers 
guidance  for  selecting  multistage  core  compressor 
design  parameters  that  have  a  high-efficiency  potential. 
Based  upon  these  results  and  some  vector  diagram  con¬ 
siderations,  figure  6  was  constructed  using  information 
from  actual  compressor  rig  tests  of  machines  having  a 
certain  degree  of  similarity  in  terms  of  tip  clearance, 
stability  requirement,  specific  flow,  and  blading.  The 
abscissa  is  the  Skoch  parameter  ( SP ),  which  is  defined  as 


SP  = 


atp 

AV'Y 


1  + 


(  i  y 

l<p-tanain> 


1/2 


where  a  is  the  solidity  (the  ratio  of  chord  to  average 
tangential  spacing);  A  is  the  blade  aspect  ratio  (the  ratio 
of  average  length  to  pitchline  chord);  and  ajn  is  the  stage 
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Figure  6.— Peak  efficiency  trends  with  various  design 
choices  for  geometric  form  and  vector  diagram. 


inlet  swirl  angle.  This  parameter,  averaged  over  all 
stages,  is  analogous  to  a  specific  speed.  The  figure  indi¬ 
cates  the  existence  of  judicious  combinations  of  geomet¬ 
ric  form  and  vector  diagram,  which  lead  to  high 
efficiencies.  This  is  discussed  by  Wisler,  Koch,  and 
Smith,  who  showed  the  competition  amongst  various  loss 
sources  as  aspect  ratio  and  solidity  are  varied.  Thus,  the 
Skoch  parameter  provides  a  basis  for  trading  geometric 
form  and  vector  diagram  (e.g.,  compactness  and  oper¬ 
ability)  to  optimize  efficiency. 

To  achieve  the  low  parts  count  and  high  efficiency 
required  of  modern  turbocompression  systems,  advanced 
technology  concepts  have  been  driven  to  high-work- 
coefficient  and  moderate-flow-coefficient  designs.  This 
trend  is  illustrated  in  figure  6  with  the  projected  advanced 
technology  goal  for  the  chosen  Skoch  parameter.  The 
advanced  technology  goal  is  a  result  of  design  choices 
orchestrated  using  the  anticipated  aerodynamic  technol¬ 
ogy  advancements.  A  discussion  of  these  design  choices 
follows. 


Four-Stage,  High-Pressure-Compressor  Aerody¬ 
namic  Technology  Demonstrator 

So  far,  no  definitive  path  has  been  identified  for  achiev¬ 
ing  the  assumed  technology  advancements.  A  research 
configuration  was  conceived  with  the  intent  of  evaluating 
how  far  the  use  of  advanced  turbomachinery  CFD  and 
three-dimensional  design  concepts  might  permit  one  to 
progress  in  realizing  high  efficiency  at  increased  loading 
levels.  The  objective  was  to  configure  a  technology 


Figure  7. — Aerodynamic  design  system. 


demonstrator  compatible  with  current  multistage  core 
compressors  for  high-bypass-ratio  turbofans.  Although 
a  specific  application  was  not  defined,  a  four-stage 
configuration  under  current  aeromechanical  constraints 
was  selected.  The  execution  of  the  aerodynamic  design 
consists  of  an  iterative  loop  between  several  computer 
programs  as  illustrated  in  figure  7.  This  methodology  is 
essentially  the  same  as  that  currently  employed  in  most 
design  offices  with  perhaps  the  exception  of  the  three- 
dimensional  inverse  design  method.  The  steps  of  the 
so-called  preliminary  design  phase  of  this  process  and 
their  interactions  are  outlined  below;  the  more  detailed 
design  phase,  including  three-dimensional  inverse  design 
and  multistage  CFD,  is  described  in  part  II. 

Meanline  performance  synthesis . — Four  parameters 
were  established  at  the  outset.  The  corrected  specific  flow 
(flow  per  unit  annulus  area)  at  the  first  rotor  inlet  was 
40.0  lbm/s-ft2  corresponding  to  an  average  axial  inlet 
Mach  number  of  0.6,  which  is  consistent  with  advanced 
core  compressors.  A  maximum  rotor  corrected  tip  speed 
of  1500  ft/s  was  assumed  based  on  turbine  stress  consid¬ 
erations  for  a  high-bypass-ratio  turbofan.  An  average 
Skoch  parameter  in  the  range  of  1 .5  to  2.0  was  selected  for 
high  efficiency  and  compactness  as  suggested  from 
figure  6.  With  a  casing  diameter  of  21 .08  in.  and  a  radius 
ratio  of  0.528  at  the  first  rotor  inlet,  the  mass  flow 
corresponding  to  the  prescribed  specific  flow  was 
70.0  lbm/s.  The  design  can  be  scaled  to  different  mass 
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TABLE  I.— COMPARISON  OF  STAGE  LOADING  AND 
POLYTROPIC  EFFICIENCY _ 


Stage 

Cumulative 

total-pressure 

ratio, 

PR„ 

State  of  the  practice 

Advanced  technology 

Effectivity 

ratio, 

E 

Polytropic 

efficiency, 

percent 

Effectivity 

ratio, 

E 

Poly  tropic 
efficiency, 

percent 

1 

2.372 

0.894 

0.890 

0.755 

0.923 

2 

4.645 

.935 

.883 

.779 

.920 

3 

8.213 

.982 

.864 

.820 

.919 

4 

12.00 

.883 

.899 

.730 

.927 

Overall 

12.00 

0.883 

0.920 

flow  rates  if  required.  Most  of  the  remaining  characteris¬ 
tics  resulted  from  the  objective  of  designing  a  multi¬ 
stage  compressor  with  a  high-efficiency  potential  at  an 
elevated  loading  level. 

The  initial  step  in  the  preliminary  design  of  the  four- 
stage,  high-pressure  compressor  (HPC)  was  taken  by 
using  a  meanline  program  incorporating  the  loss  and  peak 
pressure-rise  models  previously  presented.  This  program 
also  includes  estimates  ofhub  and  casing  vector  diagrams 
assuming  an  isentropic  simple  radial  equilibrium.  A 
conceptual  design  using  state-of-the-practice  technology 
and  current  design  practices  was  executed  first.  Then, 
assuming  the  postulated  technology  advancements,  an 
assessment  of  the  performance  potential  of  this  design 
was  made. 

The  meanline  program  provided  the  average  loss,  aero¬ 
dynamic  blockage,  and  peak  pressure-rise  capability  of 
each  stage  of  the  four-stage  configuration.  A  stagewise 
distribution  of  effectivity  was  prescribed  based  on  the 
philosophy  that  the  operating  points  of  the  front  and  back 
stages  would  primarily  pivot  about  the  two  middle  stages 
as  the  compressor  is  throttled  along  a  typical  engine 
operating  line.  In  addition,  the  averaged  effectivity  of  the 
compressor  was  selected  to  be  consistent  with  a  positive 
stability  margin  of  about  5  percent  at  the  aerodynamic 
design  speed.  The  required  stability  margin  depends  on 
the  intended  application  and  availability  of  practical 
stability-management  technologies.  Stage  effectivity 
distributions  are  given  in  table  I,  which  shows  the 
two  middle  stages  to  be  the  set  near  their  peak  loading 
capabilities. 

Of  the  several  flow-path  shapes  (constant  casing, 
constant  meanline,  and  constant  hub)  investigated,  a 
constant-diameter  casing  provided  the  best  balance 
between  high  pressure  ratio  and  good  efficiency.  Because 
of  the  need  to  limit  the  exit  rim  speed,  the  casing  diameter 
was  reduced  at  the  aft  end  of  the  compressor.  A  discharge 
Mach  number  of  about  0.35  was  assumed  based  on 


considerations  of  matching  with  an  additional  down¬ 
stream  stage.  Stage  efficiency  was  refined  by  balancing 
stage  reaction,  inlet  swirl,  solidity,  aspect  ratio,  and  axial 
velocity  diffusion.  Orchestrating  all  these  parameters 
resulted  in  an  overall  total-pressure  ratio  of  1 2: 1  at  a  rotor 
1  (R 1 )  corrected  tip  speed  of  1 477  ft/s.  The  average  aspect 
ratio  of  the  rotors  was  0.8  and  that  of  the  stators  was  1 .0. 
Aspect  ratios  for  the  rotors  were  governed  primarily  by 
aeroelastic  stability  considerations.  The  rotor  and  stator 
average  solidities  were  about  2.0.  The  predicted  stage-by- 
stage  distribution  of  polytropic  efficiency  and  com¬ 
pressor  overall  efficiency  is  given  in  table  I.  Reflected 
in  these  efficiencies  are  the  following  rotor  tip  clearance 
levels :  0.020  in .  for  R 1 , 0.0 1 8  in .  for  R2, 0.0 1 6  in.  for  R3, 
and  0.014  in.  for  R4.  An  evaluation  of  current  and  near- 
term  core  compressor  clearance  trends  supported  the 
assumed  clearance  levels.  The  first  three  stators  are 
cantilevered  from  the  casing  with  an  average  hub  clear¬ 
ance  of  0.018  in.  These  hub  clearances  and  other  clear¬ 
ances  for  anticipated  variable  stators  were  factored  into 
theefficiency  predictions.  The  overall  performance  quoted 
for  the  four-stage  HPC  extends  from  the  guide-vane 
leading  edge  to  the  compressor  discharge. 

The  flow  path  displayed  in  figure  8  was  biased  towards 
high  work  and  high  efficiency.  A  variable-stagger  inlet 
guide  vane  (IGV)  provides  part-speed  operability.  The 
guide  vane  is  intended  to  deliver  a  nominal  swirl  distribu¬ 
tion  of  9°  counterswirl  at  the  hub,  2°  counterswirl  at 
midheight,  and  0°  at  the  casing.  This  choice  was  made  to 
control  the  rotor  1  hub  reaction  and  to  reduce  the  absolute 
Mach  number  at  the  inlet  to  the  stator  1  hub.  Note  that  the 
stator  leading  edges  are  swept  forward  near  the  endwalls 
to  tailor  the  loading  in  those  regions.  There  are  32  inlet 
guide  vanes,  54  first  stators  (SI),  74  second  stators  (S2), 
132  third  stators  (S3),  and  104  fourth  stators  (S4).  The 
rotor  blade  counts  are  26  Rl’s,  56  R2’s,  76  R3’s,  and 
80  R4’s.  Preliminary  estimates  of  off-design  perfor¬ 
mance,  assuming  current  technology,  indicate  a  potential 
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Figure  8. — High-pressure  compressor  flow  path. 


need  for  variable  SI,  variable  S2,  and  interstage  bleed 
between  the  second  and  third  stages. 

The  ultimate  performance  potential  of  the  previously 
described  state-of-the-practice  design  was  evaluated  by 
incorporating  the  anticipated  advanced  aerodynamic  tech¬ 
nologies  in  the  meanline  program.  The  results  shown  in 
table  I  indicate  a  significant  increase  in  overall  efficiency 
in  addition  to  a  higher  stability  margin  for  the  same  duty. 
Thus,  the  advanced  technology  goal  displayed  in  figure  6 
was  set. 

Throughflow . — After  selecting  the  aerodynamic 
design  point  and  verifying  that  the  basic  geometric 
parameters  derived  from  the  meanline  analysis  had  the 
potential  to  be  developed  into  a  viable  configuration 
which  satisfied  the  prescribed  technical  goals,  the  span  wise 
tailoring  of  the  vector  diagrams  was  begun  by  using  an 
adaptation  of  the  inverse  throughflow  computer  program 
described  in  reference  16.  The  tailoring  involves  an 
iterative  loop  between  an  intrablade  mean-stream  surface 
analysis  and  a  blade  geometry  generator.  The  computa¬ 
tions  are  based  on  the  streamline  curvature  method  with 
the  inclusion  of  an  axisymmetric  body  force  field 
assumed  to  act  normal  to  the  surface  formed  by  the 
stacked  camber  lines  of  each  blade  row.  Custom-tailored 
airfoils  were  employed  for  all  blade  rows.  The  airfoil 
customization  technique  used  consists  of  assuming  the 
meridional  distributions  of  total  pressure  through  the 
rotors  and  of  angular  momentum  through  the  stators. 
Then,  after  correcting  for  departures  from  the  mean 
camber  line,  airfoils  are  fitted  to  the  calculated  relative 
flow  angles  and  a  two-part,  quarter-sine- wave  thickness 


distribution.  The  stream  surface  sections  are  stacked 
along  a  specified  backbone  or  stacking  axis  to  obtain 
three-dimensional  blade  surfaces. 

A  process  of  successive  refinement  guided  by  CFD 
results  determined  the  incidence  and  deviation  angles 
used  to  generate  the  airfoil  sections.  Initial  values  were 
estimated  usingthequasi-three-dimensional,  viscous  code 
RVCQ3D  (ref.  17).  RVCQ3D  was  also  used  to  direct 
the  initial  selection  (meridional  distributions  of  angular 
momentum  or  total  pressure  within  blades)  of  airfoil 
contours.  These  contours  were  updated  by  conducting 
isolated  blade-row  viscous  analysis  using  ADPAC 
(ref.  18)  with  inflow  and  outflow  boundary  conditions 
extracted  from  the  throughflow  analysis.  The  results  from 
the  ADPAC  analysis  were  also  used  to  tailor  the  stacking 
axis  and  to  modify  vector  diagrams  to  control  spanwise 
static-pressure  gradients.  Span  wise  distributions  of  blade- 
row  losses  were  also  revised  using  the  CFD  results  as  a 
guide.  Aerodynamic  blockages  were  initially  taken  from 
the  meanline  program  and  were  successively  adjusted 
based  on  CFD  results. 

Setting  the  compressor  discharge  Mach  number  and 
total-pressure  ratio  fixes  the  overall  static-pressure  rise. 
With  the  flow  path  fixed,  the  meridional  static-pressure 
gradient  is  primarily  determined  by  the  stage  effectivity 
distributions  derived  from  the  meanline  analysis.  The 
wall  static-pressure  distributions  established  with  the 
throughflow  analysis  at  the  aerodynamic  design  point 
are  displayed  in  figure  9.  To  achieve  low  loss  and  to 
ensure  aerodynamic  stability,  it  is  generally  desirable  to 
minimize  static-pressure  gradients  along  the  endwalls. 
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Figure  1 0.— Rotor  inlet  relative  Mach  number. 


Figure  9  shows  that  the  meridional  gradients  of  static 
pressure  are  approximately  constant  through  most  of  the 
rotor  blade  rows  (see  fig.  8  for  blade  locations)  and  then 
taper  off  toward  the  trailing  edge.  The  static-pressure  rise 
is  more  rapid  in  the  front  half  of  the  stator  blade  rows  and 
quickly  tapers  off  in  the  remaining  half  of  the  blade.  For 
all  blade  rows,  preferential  radial  forces  were  introduced 
through  the  tangential  leaning  of  the  mean  camber  surface 
to  control  spanwise  gradients  of  static  pressure. 

Spanwise  variations  of  the  rotor  inlet  relative  Mach 
number  are  shown  in  figure  10  for  all  four  rotors  at  the 


Figure  1 1  .—Stator  inlet  absolute  Mach  number. 


aerodynamic  design  point.  For  rotor  1 ,  the  relative  Mach 
numbers  vary  from  high  subsonic  near  the  hub  to  1 .6  near 
the  casing.  The  two  middle  rotors  (R2  and  R3)  have 
supersonic  inlet  Mach  numbers  across  the  span,  which 
makes  stage  matching  exceptionally  difficult.  The  rotor  4 
inlet  Mach  numbers  are  approximately  sonic  across  the 
span.  The  inlet  absolute  Mach  number  distributions  for 
the  stators  are  displayed  in  figure  11.  The  stator  1  hub 
Mach  number  is  slightly  above  sonic  with  all  remaining 
stators  having  subsonic  inlet  Mach  numbers. 


NASA/TP— 2002-21 1568 


11 


Percent  of  span  from  hub  Percent  of  span  from  hub 


Figure  13. — Stator  diffusion  factors. 
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Figure  14.— Stator  inlet  and  exit  absolute  flow  angles. 
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Figure  15. — Rotor  exit  cumulative  total-pressure  ratio. 
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Figure  16. — Axisymmetric  averaged  throughflow  velocity  distribution  at  aerodynamic  design 
point. 


The  span  wise  distributions  of  the  diffusion  factor  Dfactor 
for  rotors  and  stators  at  the  aerodynamic  design  point  are 
shown  in  figures  1 2  and  1 3,  respectively.  These  diffusion 
factors  are  locally  high,  reflecting  the  increased  stage 
aerodynamic  loadings.  The  maximum  diffusion  factors 
are  in  excess  of  0.65  and  occur  near  the  endwalls.  How¬ 
ever,  the  average  values  are  not  far  from  levels  for  which 
acceptable  performance  has  been  demonstrated.  The 
stator  inlet  and  exit  absolute  flow  angles  are  plotted  in 
figure  14.  On  the  average,  approximately  50°  of  turning 
was  achieved  in  the  stators.  The  intent  was  to  overturn 
near  the  endwalls  to  accommodate  secondary  flows. 

Counterswirl  (swirl  opposite  the  direction  of  rotor 
rotation)  was  prescribed  for  S2  and  S3  to  control  the 
spanwise  distributions  of  rotor  (R3  and  R4)  aerodynamic 
loading,  as  defined  by  the  diffusion  factor.  This 
counterswirl  also  increases  the  stage  reaction,  which  may 
have  a  favorable  impact  on  stage  performance.  Spanwise 
distributions  of  cumulative  rotor  exit  total-pressure  ratios 
are  displayed  in  figure  15.  These  distributions  are  com¬ 
patible  with  the  prescribed  stator  exit  swirl  profiles.  The 
intent  was  to  propagate  a  hub-strong,  total-pressure 
profile  within  the  multistage.  Such  a  profile  was  thought 
to  be  favorable  with  respect  to  aerodynamic  stability  as 
the  machine  is  throttled.  Overall,  the  vector  diagrams 
were  tailored  under  the  guidance  of  CFD  results  to  pro¬ 
vide  a  beneficial  relief  to  the  high-loss  and  high-loading 
endwall  regions. 


The  initial  blading  was  generated  and  refined  using 
isolated  blade-row  CFD  analysis.  The  axisymmetric 
averaged  throughflow  velocity  distribution  at  the  aerody¬ 
namic  design  point  is  presented  in  figure  16.  This  was 
processed  from  the  results  of  a  three-dimensional  multi¬ 
stage  simulation  employing  a  mixing  plane  approach 
within  the  ADPAC  code.  The  results  serve  to  identify  in 
the  design  potential  weak  points  that  are  in  need  of  further 
refinements.  Evident  from  the  figure  is  that  in  the  hub 
region  of  R3,  there  is  weak  flow  caused  by  the  high 
aerodynamic  loading  being  requested  and  by  aerody¬ 
namic  mismatching  between  blade  rows.  In  addition,  the 
rotor  tips  show  flow  weaknesses  associated  with  leakage 
flow  development  in  a  high-loading  region  with  shocks. 
Based  on  the  results  from  the  multistage  ADPAC  simula¬ 
tion,  the  overall  performance  potential  of  this  preliminary 
design  configuration  was  predicted  to  be  a  pressure 
ratio  of  12:1  at  a  polytropic  efficiency  of  89  percent.  To 
further  appreciate  the  level  of  technical  advancement 
being  sought,  figure  17  shows  the  fully  developed  com¬ 
pressor  goal  relative  to  established  performance  trends. 
This  design  sits  at  an  average  work  factor  of  0.44.  It  can 
be  observed  that  the  preliminary  design  results  are  in  line 
with  current  trends  and  that  achievement  of  the  goal 
would  indicate  a  substantial  advancement  over  the  state  of 
the  practice  if  not  over  the  state  of  the  art.  What  is  not 
explicitly  stated  is  the  suitability  of  such  a  design  for  a 
practical  engine  product.  Further  research  and  develop¬ 
ment  are  required  to  ascertain  the  postulated  aerodynamic 
technology  advancements. 
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Figure  17.— Core  compressor  technology  assessment  in  terms  of  work  factor  and 
polytropic  efficiency.  UHBR,  ultra-high  bypass  ratio;  EIS,  entry  into  service. 


Part  II:  Three-Dimensional  Blading 
Development 

Three-Dimensional  Inverse  Method 

The  primary  prescribed  quantities  in  the  three-dimen¬ 
sional  inverse  method  are  the  blade  stacking  line/stack  = 
/stack(r)>  t'ie  blade  thickness  distribution  T  =  7(r,z),  and 
the  blade  pressure-loading  distribution  A p  =  A p(r,z). 
Herein,  blade  pressure  loading  is  defined  as  the  pressure 
difference  between  the  blade  upper  and  lower  surfaces  at 
fixed  axial  positions.  For  a  given  set  of  inputs,  the  three- 
dimensional  inverse  method  computes  the  corresponding 
blade  mean  camber  surface  /  =  f(r,z ).  Clearly,  the 
blade  geometry  corresponding  to  prescribed  values  for 
\fs tack>  T,  A p]  is  not  guaranteed  to  have  the  optimum 
performance  or  to  be  aeromechanically  acceptable.  The 
challenge  is  to  pick  these  quantities  to  arrive  at  a  satisfac¬ 
tory  blade. 

The  three-dimensional  inverse  method  is  based  on  the 
solution  of  the  Navier-Stokes  equations  using  the  robust, 
finite-volume,  time-marching  scheme  of  Jameson, 
Schmidt,  and  Turkel  (ref.  19).  Viscous  effects  are  mod¬ 
eled  using  the  method  proposed  by  Denton  (ref.  20)  for 
turbulent  flows  as  adapted  by  Damle  (ref.  21).  During  the 
time-marching  process,  fluid  is  allowed  to  cross  the  blade 
surfaces,  and  a  pressure-jump  condition  (blade  pressure 
loading)  is  imposed  across  the  blade  surfaces.  The  “flow- 
tangency”  condition  along  the  blade  surfaces  is  then  used 
to  update  the  blade  geometry.  Endwall  clearances  are 
modeled  by  assuming  periodicity  within  the  clearance 


gap.  The  INV3D  computer  code  employed  in  this  study 
can  run  in  either  the  standard  analysis  mode  or  the  inverse 
mode. 

The  successful  strategy  employed  in  the  previously 
published  work  on  this  three-dimensional  inverse  method 
(ref.  22)  was  used  to  pick  the  prescribed  quantities 
[/stack’  that  would  improve  the  present  transonic 

rotor  blade  designs.  The  strategy  is  as  follows: 

1.  Start  with  blade  geometry  obtained  with  the  tradi¬ 
tional  design  tools  (e.g.,  meanline  and/or  throughflow 
method).  This  geometry  is  termed  the  “original  design.” 
Then  switch  to  the  three-dimensional  inverse  method  to 
improve  the  geometry,  using  the  original  design  as  the 
initial  guess. 

2.  Input  to  the  three-dimensional  inverse  method  the 
blade  tangential  thickness  distribution,  which  has  been 
kept  the  same  as  the  original  design  to  satisfy  some 
structural  constraints  and  to  reduce  the  number  of  degrees 
of  freedom. 

3.  Adjust  the  blade  pressure  loading  A p(r,z)  to  intro¬ 
duce  geometric  features  that  enable  the  control  of  local 
flow  structures  (shock,  secondary  flows,  boundary  lay¬ 
ers,  etc.)  linked  to  a  performance  index.  The  choice  of  the 
blade  pressure  loading  is  a  three-dimensional  approach. 
At  a  given  spanwise  station,  the  blade  pressure-loading 
distribution  is  used  to  manage  the  local  aerodynamics  and 
to  tailor  the  spanwise  vector  diagrams: 

£rApdA^m[(rVQ)[r(rVe)le]  (1) 
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where  the  subscripts  te  and  le  are  trailing  edge  and  leading 
edge;  Ae  is  the  tangential  projection  of  the  blade  surface 
area;  m  is  the  mass  flow  rate;  and  Vg  is  the  mass-averaged 
tangential  velocity.  For  a  rotor,  specifying  the  blade 
pressure  loading  at  every  spanwise  station  is  the  same  as 
prescribing  the  spanwise  distribution  of  total  temperature 
rise  across  the  rotor.  This  is  the  usual  two-dimensional 
thinking  when  atwo-dimensional  inverse  method  is  used. 
With  a  three-dimensional  method,  an  additional  control  is 
the  spanwise  variation  of  the  blade  pressure  loading.  As 
will  be  described  shortly,  this  is  used  to  adjust  the  orien¬ 
tation  of  the  passage  shock  as  seen  in  a  meridional  plane. 

4.  Modify  as  needed  the  prescribed  blade  stacking  line 
so  as  to  restrain  potential  endwall  corner  separations  and 
to  limit  excessive  blade  twisting. 

The  blade  design  procedure  proposed  herein  is  very 
compatible  with  the  inverse  throughflow  method  used  in 
part  I  to  generate  the  original  blades.  In  particular,  the 
design  intents,  as  described  by  the  spanwise  profiles  of 
axisy mmetric  averaged  pressure,  temperature,  and  angles, 
are  fed  into  the  inverse  method  as  inflow  and  outflow 
boundary  conditions.  One  of  the  main  advantages  of 
using  a  three-dimensional  inverse  method  is  that  the 
critical  physics  of  the  complex  aerodynamic  interactions 
are  directly  accounted  for  rather  than  being  patched  in  a 
somewhat  convoluted  fashion. 


Reblading  Using  INV3D 

The  front  two  stages,  including  the  inlet  guide  vane,  of 
the  four-stage  configuration  were  extracted  for  a  more 
refined  aerodynamic  design.  A  meridional  cross  section 
of  the  flow  path  is  shown  in  figure  18.  The  following  are 
the  aerodynamic  design  requirements  for  the  two  stage:  a 
corrected  flow  of  70  lbm/s,  an  overall  total-pressure  ratio 
of  4.645:1,  and  a  rotor  1  corrected  tip  speed  of 
1477  ft/sec.  From  part  I,  the  state-of-the-practice 
efficiency  potential  of  this  configuration  was  estimated  to 
be  89  percent  polytropic.  The  intended  axisymmetric 
profiles  and  design  considerations  for  the  original  blades 
are  found  in  part  I. 

In  this  section,  the  redesign  of  the  second  rotor  (R2)  is 
discussed  in  detail  for  three  important  reasons.  First,  the 
second  rotor  is  the  most  challenging  because  the  incom¬ 
ing  relative  Mach  number  is  supersonic  from  hub  to  tip. 
Second,  the  relative  aerodynamic  loading  based  on  the 
effectivity  parameter  E  is  the  highest.  Third,  this  rotor  is 


Figure  18— Flow  path  for  two-stage  high-pressure 
compressor. 


also  critical  in  terms  of  proper  stage  matching  because  of 
possible  strong  interactions  with  the  upstream  stator. 

The  design  intents  are  obtained  from  the  throughflow 
model,  which  is  updated  based  on  information  supplied 
by  the  multistage  CFD  simulations  to  be  described  in  the 
section  Multistage  Performance  Evaluation.  It  is  not 
possible  to  get  an  exact  match  between  the  present 
throughflow  and  the  axisymmetric  average  of  the  CFD 
solution  because  of  inconsistencies  in  formulation. 
Nevertheless,  an  acceptable  match  was  achieved.  At  the 
inlet,  the  spanwise  distributions  of  total  pressure,  total 
temperature,  and  flow  angles  are  specified.  The  design 
intent  spanwise  distribution  of  static  pressure  is  pre¬ 
scribed  downstream.  Since  the  present  formulation  of  the 
inverse  problem  called  for  the  specification  of  the  blade 
pressure  loading,  which  is  directly  related  to  the  spanwise 
distribution  of  total  temperature  rather  than  to  the  total 
pressure,  adjustments  to  the  magnitude  of  the  blade 
pressure  loading  were  sometimes  required  to  match  the 
design  intent  exit  total-pressure  distribution  and  the  mass 
flow  rate. 

Forthis  reblading  exercise,  the  meridional  envelopes  of 
the  blade  were  fixed.  The  blade  pressure  loading  was 
carefully  constructed  with  the  following  scheme.  Figure 
19  shows  the  contour  of  the  pressure  loading  of  the 
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Figure  19. — Blade  pressure  loading. 


NAS  A/TP— 2002-2 1 1 568 


15 


Figure  20. — Axial  loading 
distribution. 

original  blade  at  the  design  intent  backpressure  (left 
contour  plot).  Two  distinct  discontinuities  spanning  from 
hub  to  tip  are  clearly  indicated.  The  axial  distribution  of 
the  blade  pressure  loading  superimposed  on  the  Mach 
number  contours  in  the  blade-to-blade  plane  at  the  same 
span  location  (near  tip  region)  is  shown  in  figure  20.  Note 
that  the  two  near  discontinuities  in  the  blade  loading 
distribution  conform  to  the  location  of  the  passage  shock. 
The  front  discontinuity  corresponds  to  the  passage  shock 
impinging  on  the  pressure  surface  of  the  blade  above  it, 
whereas  the  back  discontinuity  corresponds  to  the  pas¬ 
sage  shock  impinging  on  the  suction  surface  of  the  blade 
below  it.  Looking  back  to  figure  1 9,  note  that  the  passage 
shock  is  rather  strong  at  the  upper  half  of  the  blade  and 
hence  there  exists  a  potential  for  tailoring  the  pressure¬ 
loading  distribution  to  weaken  it,  as  has  been  demon¬ 
strated  in  previous  work  (ref.  23). 

Another  important  consideration  is  the  placement  and 
orientation  of  the  passage  shock  as  seen  in  the  meridional 
plane.  At  the  design  point,  the  passage  shock  should 
reside  inside  the  blade  passage,  and  its  orientation  should 
be  such  that  it  has  maximum  obliquity  relative  to  the 
incoming  flow  to  minimize  shock  loss  (shock  sweep). 
This  suggests  that  the  passage  shock  near  the  tip  region 
should  be  placed  as  far  back  as  possible  (i.e.,  the  passage 
shock  impinging  on  the  suction  surface  near  the  trailing 
edge),  whereas  it  should  be  as  far  forward  as  possible  at 
the  other  end  (the  passage  shock  impinging  on  the  pres¬ 
sure  surface  near  the  leading  edge).  After  the  desired 
shock  position  in  the  blade  passage  is  selected,  it  is 
possible  to  estimate  where  the  passage  shock  intersects 
the  blade  surfaces  based  on  the  geometry  of  the  blade 
passage  (the  blade  stagger  angle  and  the  blade  spacing). 

There  are  a  number  of  different  aspects  to  the  selection 
of  pressure  loading  for  transonic  blades.  First,  the  selec¬ 
tion  of  loading  shape  and  magnitude  in  the  entrance 
region  of  the  blade  translates  into  tailoring  the  aerody- 


Figure  21. — Passage  shock  positioning  in  meridional 
plane. 


namic  surfaces  for  the  proper  supersonic  wave  pattern 
and  throat  area.  This  tailoring  is  a  means  for  achieving  the 
required  mass  flow  rate  with  a  started  shock  configura¬ 
tion.  Additionally,  the  selection  of  shock  position  on 
suction  and  pressure  surfaces  as  well  as  its  strength  must 
be  carefully  crafted  to  obtain  low  losses.  After  the  shock 
position  in  the  meridional  plane  has  been  selected,  one  of 
the  loading  discontinuities  is  fixed  at  that  span  position 
(points  A  and  B  in  fig.  2 1 ).  The  other  discontinuity  can  be 
placed  to  introduce  some  obliquity  of  the  shock  in  the 
blade-to-blade  plane.  As  part  of  the  selection  of  the 
meridional  shock  position,  the  region  encompassing  the 
passage  shock  must  also  be  determined.  The  tip  of  the 
blade  sets  the  first  spanwise  location  and  the  second  is 
chosen  to  be  where  the  inlet  relative  Mach  number  is 
slightly  above  unity.  This  positioning  divides  the  blade 
into  the  low  supersonic-subsonic  region  where  smooth 
loading  shapes  can  be  used  without  producing 
discontinuities  in  the  blade  and  the  supersonic  flow 
portion  where  the  discontinuities  in  loading  must  be 
prescribed  to  obtain  “smooth”  blade  shapes  (fig.  21). 

The  strength  of  the  shock  can  be  controlled  somewhat 
through  the  specification  of  the  severity  of  the  pressure¬ 
loading  gradient.  A  steep  gradient  will  produce  a  strong 
shock  with  a  greater  shock  loss  and  a  potentially  greater 
risk  of  shock-induced,  boundary-layer  separation.  How¬ 
ever,  there  are  limits  as  to  how  much  the  shock  can  be 
weakened  and  how  much  obliquity  can  be  introduced  in 
the  blade-to-blade  plane.  If  the  shock  is  weakened  past  a 
certain  point,  the  resulting  blade  will  have  a  “kink”  that 
may  lead  to  poor  performance  at  off-design  conditions. 
The  second  limit,  blade-to-blade  shock  obliquity,  is 
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Figure  22. — Blade-surface  pressure  and  loading  distributions  for  several  span  sections,  (a)  Near  the  hub.  (b)  Below 
50-percent  span,  (c)  Above  50-percent  span,  (d)  Near  the  tip. 


observed  when  the  two  discontinuities  have  been  placed 
too  far  from  one  another  (dashed  line  in  fig.  21).  In  such 
a  case,  the  shock  wi  11  occur  at  the  specified  location  on  the 
pressure  surface  (front  discontinuity,  see  fig.  21)  but  will 
impinge  on  the  suction  surface  at  a  location  upstream  of 
the  “specified”  back  discontinuity,  resulting  in  too  much 
loading  in  the  subsonic  diffusion  portion  of  the  blade  and 
therefore  overcambering  the  blade  in  that  region.  This 
“manual”  adjustment  of  the  blade  pressure  loading  is  a 
delicate  process  and  does  require  some  experience  to 
achieve  success. 

Figure  1 9  (right  contour  plot)  shows  the  blade  pressure 
loading  of  the  redesigned  rotor  2.  Compared  with  the 
original  design  (left  contour  plot),  the  figure  clearly 


illustrates  the  weakening  of  the  passage  shock  and  an 
increase  in  its  obliquity.  Figure  22  compares  the  detailed 
blade  pressure  loading  (bottom)  and  pressure  distribu¬ 
tions  on  the  blade  surfaces  (top)  of  the  original  design  and 
the  new  design  at  several  spanwise  stations.  Two  impor¬ 
tant  differences  between  the  two  designs  are  worth  point¬ 
ing  out.  First,  the  passage  shock  is  weakened  at  all 
spanwise  stations.  Second,  the  passage  shock  moves  a 
little  forward  below  the  50-percent-span  location  while 
being  displaced  farther  aft  above  the  50-percent-span 
location.  Recall  that  the  reorientation  of  the  passage 
shock  was  purposely  created  in  the  specification  of  the 
blade  pressure  loading  to  increase  its  obliquity  relative  to 
the  meridional  flow. 
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Figure  23.— Relative  Mach  number  contours  at  blade 
trailing-edge  plane. 


Adiabatic  efficiency,  v\tt 


Figure  24. — Comparison  of  spanwise  distribution  of 
adiabatic  efficiency  at  design  point. 

Figure  23  shows  the  flow  field  in  a  crossflow  (r,0)  plane 
behind  the  blade  trailing  edge.  The  tip  clearance-to-chord 
ratio  is  approximately  1  percent.  For  the  new  design,  there 
is  a  clear  indication  of  a  cleaner  endwall  flow  that  is 
attributed  to  the  weakened  passage  shock  and  possibly  a 
modulation  of  secondary  and  leakage  flows.  Note  that  the 
new  blade  geometry  is  twisted  differently  near  the  casing 
endwall,  whereas  the  original  blade  was  stacked  so  that 
its  leading  edge  is  swept  forward  relative  to  the  flow.  This 
new  twist  introduces  a  component  of  blade  force  that 
restricts  unfavorable  spanwise  flow  migrations.  The  fig¬ 
ure  also  shows  that  the  blade  wake  for  the  new  design  is 


Figure  25. — Comparison  of  pressure  loading. 


Mach  number 

Figure  26— Cascade  flow  near  casing. 


slightly  thinner  over  the  entire  span.  The  Navier-Stokes 
solver  ADPAC  (ref.  1 8)  was  used  to  independently  con¬ 
firm  the  improved  flow  field  and  performance  of  the  new 
design.  The  calculation  used  a  mesh  size  of 90  X48  X40  (a 
total  of  172  800  cells)  that  included  4  cells  in  the  tip 
clearance  region.  At  the  design  intent  backpressure,  the 
ADPAC  solution  predicted  an  overall  adiabatic  effi¬ 
ciency  improvement  of  2  percent  over  the  original  design 
running  at  the  same  backpressure  with  approximately  the 
same  flow  rate.  Figure  24  shows  the  spanwise  distribution 
of  adiabatic  efficiency  as  predicted  by  ADPAC.  It  can  be 
observed  that  efficiency  improvements  are  found  over  the 
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Figure  27.— Visualization  of  rotor  blades  generated  with  INV3D. 


entire  span,  with  about  1 .5  percent  over  the  core  flow  and 
over  10  percent  in  the  casing  endwall. 

The  first  rotor  (Rl)  was  also  revised  using  the  same 
pressure-loading  scheme.  Figure  25  compares  the  blade 
pressure  loading  of  the  original  design  (left  contour 
plot)  with  that  of  the  new  design  (right  contour  plot),  and 
figure  26  presents  the  blade-to-blade  flow  field  near  the 
tip.  As  in  the  second  rotor  (R2),  the  passage  shock  was 
weakened  and  slightly  reoriented  to  increase  shock  obliq¬ 
uity.  The  performance  of  Rl  (as  predicted  by  ADPAC) 
showed  an  improvement  of  1 .5  percent  in  adiabatic 
efficiency  over  the  original  design. 

In  summary,  the  transonic  rotor  blades  generated  by 
INV3D  using  the  present  pressure-loading  strategy  indi¬ 
cated  a  relatively  higher  efficiency  potential  than  the 
original  blades.  Stacking  and  thickness  distributions  were 
fixed  according  to  the  original  blades.  Perhaps  a  combi¬ 
nation  inverse  direct  method  can  be  used  to  explore  other 
degrees  of  freedom.  Preliminary  mechanical  and  struc¬ 
tural  assessments  of  the  inverse-designed  blades  thus  far 
appear  promising.  The  blades  are  curvier  than  usual, 
having  forward-swept  leading  edges  (fig.  27).  Implica¬ 
tions  for  this  design  freedom,  from  the  perspective  of 
product  development  and  total  cost,  are  beyond  our 
scope.  Similar  to  the  rotors,  a  pressure-loading  strategy 
can  be  developed  for  subsonic  stators,  but  this  is  reserved 
for  future  work.  The  stators  were  designed  as  described 
in  part  I  and  were  matched  to  the  rotors  solely  based  on 
CFD  simulations.  Additional  geometric  information  for 
all  blade  rows  is  given  in  appendix  A. 


Multistage  Performance  Evaluation 

Procedure. — To  verify  the  improved  performance 
potential  of  the  rebladed  two  stage  and  to  facilitate  stage 
matching,  the  multistage  CFD  code  APNASA,  based  on 
the  formulation  of  Adamczyk  (ref.  5),  was  employed.  The 
APNASA  code  computes  the  time-mean  flow  as  seen  by 
an  average  passage  of  a  blade  row  imbedded  in  a  multi¬ 
stage  environment.  Details  of  the  numerics,  turbulence 
model,  and  deterministic  stress  closure  are  given  by 
Adamczyk  et  al.  (ref.  24).  This  code  has  undergone 
extensive  validation  and/or  calibration  (refs.  25  and  26) 
and  has  been  found  to  be  a  reliable  predictive  tool  for 
multistage  aerodynamics.  In  executing  the  design, 
APNASA,  ADPAC,  and  INV3D  were  integrated  in  the 
iterative  design  loop.  The  procedure  is  as  follows: 

1 .  From  the  throughflow  model,  extract  the  inflow  and 
outflow  boundary  conditions  for  INV3D.  Rotor  blades 
are  designed  using  INV3D  with  the  pressure-loading 
strategy  previously  described;  the  stators  are  designed 
according  to  the  inverse  throughflow  method  outlined  in 
parti. 

2.  ADPAC  is  used  to  fine-tune  the  individual  blade-row 
designs  while  holding  the  inflow  conditions  extracted 
from  the  throughflow  model.  The  exit  profiles  from  the 
individual  blades  are  checked  against  the  design  intent 
and  the  blades  are  redesigned  as  in  step  1  should  the  exit 
profiles  or  losses  be  far  too  different  from  those  of  the 
throughflow. 
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TABLE  IL— PREDICTED  PERFORMANCE  POTENTIAL 
AT  AERODYNAMIC  DESIGN  POINT 


Inlet  correct  flow,  Wc ,  lbni/s . 71.2 

Overall  total -pressure  ratio,  PRU . 4.65:1 

Corrected  mechanical  speed,  Q ,  rpm . 16  060 

Overall  efficiency,  percent 

Adiabatic,  rj„ . 90.36 

Polytropic,  T|;> . 92.19 


3.  The  multistage  aerodynamics  of  all  the  blade  rows 
are  predicted  with  APNASA.  An  interrogation  of  the 
APNASA  solution  is  made  to  arrive  at  appropriate  revi¬ 
sions  for  the  throughflow  (i.e.,  aerodynamic  blockage, 
loss  and  turning  distributions)  model  and  the  blade  design 
strategies  (incidence,  trailing-edge  loading  distribution, 
and  stacking). 

The  above  steps  usually  converge  after  two  to  three 
design  iterations.  Note  that  step  3,  involving  mainly 
APNASA,  is  more  like  a  development  stage  rather  than  a 
design  stage.  The  throughflow  model  is  still  at  the  center 
of  the  overall  process  and  allows  all  the  blade  rows  to  be 
designed  simultaneously.  Additional  work  is  needed  to 
streamline  this  procedure  and  make  it  more  direct. 
Perhaps  a  totally  different  process  with  a  well-founded 
formalism  is  warranted.  In  addition,  the  entire  process  or 
parts  of  it  may  be  very  amenable  to  various  optimization 
schemes  that  can  directly  accommodate  aeromechanical 
and  manufacturing  constraints. 

APNASA  results . — A  sheared  H-mesh  was  constructed 
for  the  two-stage  configuration  that  included  the  inlet 
guide  vane.  The  mesh  consisted  of  61  streamwise  points 
on  each  of  the  5  blade  rows,  61  spanwise  points,  and  51 
blade-to-blade  points  for  a  total  dimension  of 626  X  6 1  X  5 1 . 
There  are  four  cells  in  the  rotor  tip  clearances  and  stator 
hub  clearances.  The  clearance  gaps  were  modeled  with  a 
periodic  boundary  condition  rather  than  with  direct 
discretization.  Different  mesh  sizes  were  investigated  to 
ascertain  the  sensitivity  of  the  results  and  to  determine  the 
adequacy  of  the  selected  mesh  size.  Simulations  were 
conducted  at  several  operating  speeds  and  backpressures. 
The  results  for  the  simulation  closest  to  the  aerodynamic 
design  point  will  be  presented  followed  by  the  overall 
characteristics  at  other  operating  conditions. 

The  predicted  performance  potential  at  the  aerody¬ 
namic  design  point  is  presented  in  table  II.  The  perfor¬ 
mance  is  determined  from  computed  circumferentially 
mass-averaged  total  pressures  and  temperatures  at  recti¬ 
fying  planes  situated  midway  between  the  individual 
blade  rows.  The  overall  performance  quoted  in  the  table 


TABLE  III.— DESIGN  POINT  STAGE 
PERFORMANCE  POTENTIAL 
PREDICTED  BY  APNASA 


Stage 

Rotor 

Stator 

Total-pressure 

Polytropic 

Total 

ratio, 

efficiency, 

pressure 

PR„ 

percent 

loss, 

A PJP„ 
percent 

1 

2.423 

95.100 

1.636 

2 

1.995 

94.420 

1.725 

Efficiency,  percent 

Figure  28. — Predicted  spanwise  variations  in  adiabatic 
efficiency  at  aerodynamic  design  point,  (a)  Overall, 
(b)  Rotor  1 .  (c)  Rotor  2. 
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Figure  29.— Relative  Mach  number  distributions  at  design  point  for  rotor  1.  (a)  9-percent 
span,  (b)  50-percent  span,  (c)  95-percent  span. 


extends  from  the  IGV  inlet  to  one  meridional  chord 
downstream  of  S2.  The  total  pressure  loss  across  the 
IGV  is  A Pt/P(  =  0.476  percent.  A  summary  of  the  stage- 
by-stage  performance  potential  is  given  in  table  III. 

At  the  design  corrected  speed  and  pressure  ratio,  the 
predicted  mass  flow  rate  is  approximately  1.7  percent 
higher  than  the  design  intent.  The  polytropic  efficiency 
significantly  exceeds  the  state-of-the-practice  value  indi¬ 
cated  in  part  I  and  is,  in  fact,  at  a  level  consistent  with 
the  four-stage  advanced  technology  goal.  The  design- 
point  overall  performance  potential  for  the  rotors  is 
exceptionally  good  in  terms  of  efficiency.  The  spanwise 
variations  of  adiabatic  efficiency  for  Rl,  R2,  and  the 
overall  two  stages  are  plotted  in  figure  28.  It  is  evident 
that  the  efficiency  of  Rl  is  relatively  high  over  the  lower 
80-percent  span  and  rapidly  drops  off  in  the  remaining 
10-percent  span  near  the  casing.  A  similar  trend  is 
observed  for  R2,  except  that  the  efficiency  levels  are 
somewhat  lower  between  20-  to  50-percent  span.  The 
distribution  of  efficiency  for  the  overall  two  stages,  ex¬ 
tracted  one  meridional  chord  aft  of  S2,  indicates  a  good 
distribution  across  most  of  the  span  with  a  gradual  dropoff 
near  the  casing.  Additional  spanwise  profiles  of 
axisymmetric  averaged  flow  properties  are  given  in 
appendix  B. 

The  APNASA  results  at  the  aerodynamic  design  point 
were  interrogated  to  identify  the  major  flow  features 
responsible  for  the  aerodynamic  character  of  the  blading 
and  to  synthesize  their  relative  impact  on  performance. 
The  cascade  plane  relative  Mach  number  distributions  for 
Rl  at  near  hub  (9-percent  span),  midspan,  and  near  tip 
(95-percent  span)  are  shown  in  figure  29.  Near  the  hub. 


Pressure  surface 


Figure  30. — Entropy  distribution  at  trailing-edge  plane  of 
rotor  1 . 

a  rapid  diffusion  from  about  40-percent  chord  without 
any  shock  is  observed  within  the  passage,  which  is 
consistent  with  the  intended  pressure  loading  shown  in 
figure  25  and  does  not  appear  to  adversely  impact  the 
viscous  layers  because  of  the  favorable  spanwise  sweep¬ 
ing  of  low-momentum  fluid  near  the  hub.  At  midspan, 
an  oblique  shock  exists  in  the  passage  at  about  40-percent 
chord,  and  the  subsonic  diffusion  in  the  aft  section  thick¬ 
ens  the  suction-surface  boundary  layer  near  the  trailing 
edge.  The  low-momentum  fluid  in  this  region  centrifuges 
outward  and  accumulates  near  the  tip.  A  cesspool,  fed 
by  this  process  and  augmented  by  the  leakage  flow,  is 
observed  at  the  blade  tip  near  the  pressure  side  comer. 
Further  evidence  of  this  cesspool  is  shown  in  figure  30, 
which  corresponds  to  the  entropy  distribution  at  the 
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Figure  31 . — Absolute  Mach  number  distribution  at 
trailing-edge  plane  of  stator  1 . 


trailing-edge  plane  of  Rl.  Near  the  tip,  note  the  existence 
of  a  double  shock  (evident  in  fig.  29)  that  is  inconsistent 
with  the  intent  shown  in  figure  25.  This  double  shock  is 
perhaps  an  indication  of  a  deficiency  in  the  viscous 
treatment  of  the  inverse  method  and  is  under  investiga¬ 
tion.  Nevertheless,  the  flow  structure  indicates  a  rela¬ 
tively  clean  flow  conducive  toahigh-efficiency  potential. 

The  absolute  Mach  number  distribution  at  the  trailing 
edge  of  S 1  is  shown  in  figure  3 1 .  Note  the  low-momentum 


region  at  the  hub  resulting  from  the  interaction  of  the  hub 
leakage  with  the  endwall  flow.  There  is  indication  of  a 
thin  wake  in  the  core  region.  Because  of  the  hub  clearance 
and  the  local  orientation  (e.g.,  sweep  and  dihedral)  of  the 
blading  near  the  endwalls,  there  is  a  tendency  for  low- 
momentum  fluid  on  the  blade  surfaces  to  migrate  into  the 
core  region,  hence  the  appearance  of  a  relatively  thinner 
wake  near  the  endwalls. 

Figure  32  shows  the  distribution  of  the  relative  total 
pressure  at  several  cross-passage  planes  from  the  leading 
edge  to  the  trailing  edge  of  R2.  Evident  is  the  develop¬ 
ment  of  the  boundary  layer  on  the  suction  surface  along 
with  the  evolution  of  a  low-energy  cesspool  associated 
with  shock-endwall  viscous  interaction  and  tip  clearance 
leakage  in  the  casing  endwall  region.  A  band  of  low 
relative  total  pressure  can  be  observed  at  the  hub,  indicat¬ 
ing  a  weak  flow  due  to  excessive  diffusion  imposed  by 
the  blading.  Although  not  explicitly  shown  here,  a  strong 
interaction  exists  between  the  bow  shock  emanating  at  the 
leading  edge  of  R2  and  the  suction-surface  boundary 
layer  at  the  trailing  edge  of  S 1 .  This  interaction  leads  to  a 
local  undertuming  of  the  stator  core  flow.  However,  there 
does  not  seem  to  be  any  adverse  effect  on  the  swallowing 
capacity  of  R2  or  on  the  overall  matching  of  the  two 
stages. 


Accumulation  of 
low-momentum  fluid 

Losses  from 
oblique  shock 


Trailing  edge 

^-Trailing-edge 
/  J  wake 


Leading 

edge 


Figure  32.— Relative  total-pressure  distribution  on  several  cross¬ 
passage  planes  for  rotor  2  (courtesy  of  ASE  Technologies,  ref.  27). 
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Figure  33. — Absolute  total-pressure  distribution  on 
cross-passage  planes  through  stator  2  (courtesy 
of  ASE  Technologies,  ref.  27). 


The  distribution  of  absolute  total  pressure  at  several 
cross-passage  planes  from  the  leading  edge  to  the  trailing 
edge  of  S2  is  presented  in  figure  33.  Clearly  evident  is  the 
development  of  the  hub  leakage  flow  and  its  interaction 
with  the  endwall  blading.  A  hub-strong,  total-pressure 
profile  exists  at  the  inlet  to  S2  and  is  rapidly  degraded  by 
this  process. 

A  loss  audit  was  performed  on  R1  to  identify  sources  of 
entropy  generation  and  their  relative  contribution  to  the 
rotor  overall  inefficiency.  The  approach  selected  for 
the  loss  audit  was  to  create  the  postprocessing  tool 
DISECT-N  (ref.  27)  that  can  decompose  the  spatial 
domain  of  a  blade  passage  into  zones  encompassing 
selected  flow  features  and  entropy  generation  mecha¬ 
nisms.  Entropy  fluxes  are  summed  over  the  boundaries  of 
each  zone  and  the  results  used  as  a  direct  measure  of  the 
entropy  production  within  that  zone.  Note  that  this  tacitly 
assumes  a  negligible  thermal  entropy  production.  Nine 
zones  (figs.  34  and  35)  were  selected  based  on  interroga¬ 
tions  of  the  APNASA  results;  the  targeted  entropy 
generation  mechanisms  are  cursorily  described  as 


Figure  34.— Zonal  decomposition  of  rotor  1 
flow  domain  for  assessment  of  entropy 
production. 


Figure  35. — Additional  zones  for 
assessment  of  rotor  1  entropy 
production. 


Zone  01 :  Casing  endwall  region — focus  on  clearance 
leakage  and  shock-endwall  boundary-layer  interactions 
Zone  02:  Casing  endwall  region,  low-momentum 
region,  including  shock  interaction  effects — focus  on 
pressure  surface 

Zone  03:  Hub  endwall  region,  low-momentum 
region — focus  on  pressure  surface 
Zone  04:  Suction-surface,  low-momentum  region 
encompassing  passage  shock-boundary-layer  inter¬ 
action — focus  on  endwall  region 
Zone  05:  Suction-surface,  low-momentum  region 
encompassing  shock-boundary-layer  interaction — focus 
on  core  flow  region 

Zone  06:  Downstream  region,  mixing  losses 
Zone  07:  Outer  region,  passage  shock — focus  on 
suction  side 
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2  4  6  8 


Zone 

Figure  36.— Relative  distribution  of  entropy  production 
for  rotor  1 . 


Zone  08:  Outer  region,  passage  shock — focus  on 
pressure  side 

Zone  09:  Inner  region,  passage  shock — focus  on  core 

The  loss  audit  results  are  given  in  figure  36,  where  the 
entropy  production  from  each  zone  is  indicated  in  terms 
of  its  relative  contribution  to  the  overall  entropy  produc¬ 
tion.  Also  shown  for  comparison  is  a  similar  loss  audit  for 
the  initial  R1  designed  using  an  inverse  throughflow 
method.  As  expected,  the  end  wall  zone  01  has  the  largest 
loss  contribution  with  approximately  30  percent  of  the 
total.  Zone  02  accounts  for  a  significant  20  percent  of  the 
total  losses  by  capturing  the  transport  of  low-momentum 
fluid  and  its  interaction  with  the  passage  shock  near  the 
pressure  surface.  The  processes  within  these  two  endwall 
zones  account  for  a  little  over  50  percent  of  the  overall 
losses.  A  significant  44  percent  of  the  total  loss  is  attrib¬ 
uted  to  the  suction-surface  viscous  dissipation  within 
zones  04  and  05  and  the  combined  outboard  shock  loss 
from  zones  07  and  08.  The  downstream  mixing  losses  are 
on  the  order  of  1  percent  of  the  total  but  do  not  cover  full 
mixing  of  the  downstream  flow  field.  Note  that  relative  to 
the  inverse  throughflow  design,  the  prescribed  pressure¬ 
loading  design  has  significantly  lower  shock  and  shock¬ 
boundary-layer-related  losses. 

Further  information  about  the  entropy  production 
mechanisms  in  the  casing  endwall  region  (zone  01, 
15-percent  immersion  from  the  casing)  was  gained  by 
subdividing  this  region  into  seven  smaller  zones  encom- 


0.825  -  Percent  l  go  l  90  \  )  100 

:  speed  /20/4?/ /10/2.3/ 95  /0/0/ 

:  /IGV  reset/  |  '  /5/1.2/ 

0.800  -  SI  reset  70 

I  /30/7/ 

Symbols  denote  APNASA 
0.775  -  predictions  at  nominal  IGV 

and  SI 


20  30  40  50  60  70  80 

Corrected  mass  flow,  m,  Ibm/s 

Figure  37. — Predicted  overall  performance  characteristics, 
(a)  Total-pressure  ratio,  (b)  Adiabatic  efficiency. 


passing  flow  features  thought  to  be  critical  and  then 
performing  a  loss  audit  for  the  region.  This  audit  indicated 
that  tip  clearance  leakage  and  viscous  dissipation  along 
the  endwall  suction-surface  comer  contribute  to  approxi¬ 
mately  40  percent  of  the  total  casing  endwall  entropy 
production.  The  passage  shock  and  its  interaction  with  the 
viscous  endwall  flow  under  the  influence  of  the  leakage 
vortex  accounts  for  about  30  percent  of  the  endwall 
entropy  production.  Although  the  spanwise  transport  of 
low-momentum  fluid  feeding  the  cesspool  near  the  pres¬ 
sure  surface  comer  directly  contributes  a  relatively  small 
amount  to  the  endwall  entropy  production,  its  indirect 
impact  on  performance  is  linked  to  the  downstream 
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mixing  loss  and  the  increased  aerodynamic  blockage. 
Thus,  with  this  audit,  various  redesign  strategies  can  be 
conceived  and  tested  via  numerical  simulation  to  assess 
their  relative  benefits. 

APNASA  was  used  to  simulate  the  off-design  perfor¬ 
mance  of  the  two-stage  configuration  at  100-,  95-,  and 
90-percent  rotational  speeds  with  nominal  IGV  and  vane 
settings.  The  geometry  was  frozen  as  designed  without 
accounting  for  structural  deformation  effects  with  vary¬ 
ing  operating  conditions.  The  results  of  the  simulation 
were  used  to  construct  stage  characteri  sties  for  each  of  the 
two  stages  and  these  were  subsequently  stacked  to  predict 
the  overall  performance  characteristics  with  various 
IGV  and  vane  resets.  The  overall  adiabatic  efficiency 
predicted  by  APNASA  was  degraded  by  1 .46  percent, 
based  on  engineering  judgment  reflective  of  hardware 
quality  shortfalls  and  other  “X-factors”  that  may  compro¬ 
mise  the  CFD  predictive  capability.  Figure  37  shows  the 
overall  characteristics  (total-pressure  ratio  and  adiabatic 
efficiency)  derived  from  the  stage-stacking  procedure  at 
the  indicated  resets  and  compares  them  with  APNASA 
predictions.  Since  there  is  no  definitive  application  for 
this  two-stage  group,  the  adequacy  of  the  characteristic 
map  cannot  be  absolutely  assessed.  Nevertheless,  a 
potentially  adequate  operating  range  seems  to  be  avail¬ 
able  for  further  aerodynamic  development. 


Conclusions 

In  part  I,  very  aggressive  performance  goals  were 
defined  by  a  rational  perturbation  of  basic  aerodynamic 
technology  elements  of  current  design  systems.  A  poten¬ 
tial  for  a  substantial  expansion  of  the  aerodynamic  design 
and  operability  space  was  shown.  The  research  at  NASA 


Glenn  on  practical  means  to  effect  these  aerodynamic 
advancements  is,  in  reality,  still  in  its  infancy.  The  chal¬ 
lenge  of  discovering  and  developing  new  technology  in 
a  multistage  turbocompressor  environment  requires  a 
well-conceived  plan  of  attack.  Nevertheless,  achieve¬ 
ment  of  these  advancements  would  give  the  design  engi¬ 
neer  the  freedom  to  innovate. 

A  preliminary  design  of  a  four-stage  compressor  con¬ 
figuration  was  established  for  advanced  technology 
development.  This  design  was  biased  toward  high  effi¬ 
ciency  at  high  loading.  Three-dimensional  blading  and 
spanwise  tailoring  of  vector  diagrams  were  employed 
under  the  guidance  of  computational  fluid  dynamics 
(CFD)  to  control  the  aerodynamics  of  the  high-loaded 
endwall  regions.  The  preliminary  results  were  encourag¬ 
ing  and  the  front  two  stages  were  extracted  for  further 
aerodynamic  trimming  using  athree-dimensional  inverse 
design  method. 

In  part  II,  the  design  principles  for  constructing  blade- 
surface,  pressure-loading  distributions  to  be  applied  in  a 
three-dimensional  inverse  design  method  were  presented. 
This  method  produced  blading  having  a  high-efficiency 
potential  relative  to  current  practice.  In  addition,  the 
inverse-designed  blades  were  successfully  adapted  to  the 
multistage  environment  using  the  APNASA  code.  An 
overall  efficiency  potential  farexceeding  current  practice 
was  demonstrated  with  the  APNASA  CFD. 

In  spite  of  the  significantly  higher  aerodynamic  load¬ 
ings,  advanced  CFD-based  tools  were  able  to  effectively 
guide  the  design  of  a  very  efficient  two-stage  compressor 
under  state-of-the-art  aeromechanical  constraints. 
Adequate  operability  for  aerodynamic  research  and  devel¬ 
opment  was  predicted.  This  two-stage  configuration  can 
be  used  for  further  aerodynamic  technology  research  and 
development  to  meet  the  goals  established  in  part  I. 
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Appendix  A 
Geometry 


TABLE  IV.— FLOW-PATH  COORDINATES:  DESIGN  AERODYNAMIC  INTENT 


Hub 

Casing  1 

Hub 

Casing 

r, 

z. 

r, 

z. 

r, 

z. 

r, 

in. 

in. 

in. 

in. 

in. 

in. 

in. 

in. 

-10.000 

4.993 

-10.00 

10.540 

8.439 

8.958 

8.071 

10.540 

-9.000 

4.993 

-9.000 

10.540 

8.628 

8.987 

8.347 

10.540 

-8.000 

4.993 

-8.000 

10.540 

8.919 

9.021 

8.695 

10.540 

-7.044 

4.993 

-7.017 

10.540 

9.211 

9.039 

9.043 

10.540 

-6.090 

4.993 

-6.040 

10.540 

9.504 

9.051 

9.392 

10.540 

-5.136 

4.993 

-5.082 

10.540 

9.796 

9.064 

9.740 

10.540 

-4.181 

4.993 

-4.169 

10.540 

10.088 

9.085 

10.088 

10.540 

-3.225 

4.993 

-3.262 

10.540 

10.224 

9.099 

10.250 

10.521 

-2.717 

5.015 

-2.766 

10.540 

10.356 

9.116 

10.408 

10.497 

-2.417 

5.039 

-2.400 

10.540 

10.546 

9.145 

10.579 

10.470 

-2.118 

5.071 

-2.031 

10.540 

10.736 

9.177 

10.749 

10.439 

-1.820 

5.112 

-1.660 

10.540 

10.927 

9.206 

10.919 

10.406 

-1.523 

5.162 

-1.287 

10.540 

11.118 

9.224 

11.089 

10.371 

-1.228 

5.221 

-.913 

10.540 

11.311 

9.221 

11.258 

10.334 

-.605 

5.377 

-.412 

10.540 

11.604 

9.182 

11.604 

10.253 

.000 

5.569 

.000 

10.540 

11.815 

9.133 

11.813 

10.202 

.734 

5.886 

.518 

10.540 

12.023 

9.071 

12.021 

10.149 

1.420 

6.294 

1.039 

10.540 

12.229 

9.004 

12.230 

10.096 

2.094 

6.721 

1.559 

10.540 

12.437 

8.940 

12.439 

10.043 

2.777 

7.134 

2.080 

10.540 

12.647 

8.886 

12.647 

9.991 

3.489 

7.495 

2.601 

10.540 

12.842 

8.850 

12.887 

9.932 

3.869 

7.645 

3.072 

10.540 

12.981 

8.834 

13.008 

9.903 

4.258 

7.766 

3.551 

10.540 

13.120 

8.826 

13.129 

9.875 

4.677 

7.873 

4.052 

10.540 

13.259 

8.824 

13.251 

9.848 

5.100 

7.964 

4.553 

10.540 

13.399 

8.822 

13.372 

9.822 

5.524 

8.052 

5.054 

10.540 

13.538 

8.819 

13.494 

9.797 

5.947 

8.146 

5.555 

10.540 

13.785 

8.806 

13.786 

9.741 

6.367 

8.252 

6.057 

10.540 

13.979 

8.790 

13.978 

9.707 

6.561 

8.310 

6.359 

10.540 

14.172 

8.771 

14.170 

9.678 

6.754 

8.374 

6.662 

10.540 

14.365 

8.749 

14.363 

9.651 

7.052 

8.485 

6.889 

10.540 

14.557 

8.727 

14.557 

9.629 

7.347 

8.603 

7.116 

10.540 

14.751 

8.713 

14.751 

9.610 

7.642 

8.719 

7.342 

10.540 

14.962 

8.710 

14.962 

9.594 

7.940 

8.828 

7.569 

10.540 

15.466 

8.709 

15.466 

9.575 

8.245 

8.918 

7.796 

10.540 

15.970 

8.709 

15.970 

9.575 
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Figure  38. — Inlet  guide  vane,  aft  ^ 

looking  forward.  Aspect  ratio,  Figure  40.— Stator  1 .  Aspect  ratio,  1 .1 28;  number  of 

3.1 50;  number  of  blades,  32.  blades,  54. 


Figure  39. — Rotor  1 .  Aspect  ratio,  0.837;  number  of  Figure  41 .  Rotor  2.  Aspect  ratio,  0.888;  number 

blades,  26.  of  blades,  56. 
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Figure  42.— Stator  2.  Aspect  ratio,  0.923;  number  of 
blades,  74. 


Figure  43. — Four-stage  flow  path  with  first  two  stages 
extracted  for  further  aerodynamic  development. 
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Fraction  of  span  Fraction  of  span  Fraction  of  span 


Figure  47. — Rectifying  plane  of  rotor  2  exit  and  stator  2  inlet. 
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